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1 Summary 

This paper defines a model that could be used to evaluate the ef- 
fectiveness of control measures during the West Africa Ebola epi- 
demic. The introduction gives a general background of the out- 
break, an overview of the disease in context, and common modes of 
regional travel. 

A brief introduction of compartmental models is given, em- 
phasizing parameters useful in modeling. Two major factors are 
addressed: the declining effectiveness of traditional control and the 
possible benefit of increased monitoring and intermittent house- 
level quarantine. 



2 Purpose 



This model has been formulated to aid implementation of the Sierra 
Leone Government's control measures, specifically the develop- 
ment of a national surveillance program combined with rolling 
lockdown. 

Although a large scale Ebola epidemic has never been controlled 
in this manner, the hypothetical effect may be inferred from the 
fact that the host species and humans have historically cohabited 
without a widespread epidemic, or reported outbreak, in West 
Africa. 1 

Further, experience has shown that calls for international assis- 
tance are not necessarily heeded in a timely or effective manner, 
and that the capacity to stop the spread of the disease through con- 
ventional measures is not scaling with the course of the epidemic. 



1 Derek J. Taylor, Matthew J. Ballinger, 
Jack J. Zhan, Laura E. Hanzly, and 
Jeremy A. Bruenn. Evidence that 
ebolaviruses and cuevaviruses have 
been diverging from marburgviruses 
since the miocene. Peer], 2:e556, 
9 2014. ISSN 2167-8359. doi: 
10.7717/peerj.556. URLhttp: 
//dx. doi.org/10.7717/peerj .556 
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3 Background 

The first identified case died on December 6 th , 2013, in Meliandou 
Village, Gueckedou Prefecture, Guinea. 2 Before being identified 2 Baize et al., 2014 

as Zaire Ebola virus (ZEBOV), the outbreak spread to Macenta pre- 
fecture around March 7 th . The disease was identified by Medecins 
sans Frontieres (MSF) from samples taken in mid-March. In early 
April, the outbreak spread to Conakry, a major urban center and 
the capital of Guinea. In early June, outbreaks were reported in 
Telimele and Boffa prefectures in Guinea, Kailahun and Kenema 
Districts in Sierra Leone, and Lofa and Montserrado counties in 
Liberia. 

The outbreak is the largest cluster of Ebola cases on record. Al- 
though previous outbreaks of ZEBOV have seen fatality rates as 
high as 90%, a rough figure of the fatality rate of the this strain of 
ZEBOV is 60%. The incubation period is approximately five days, 
but can range from 3-21 days. The infective period is approximately 
five days. The greatest risk of infection is from very sick individuals 
and the dead. 



3.1 Ebola in West Africa 



Possible reservoirs include fruit bats of the species Hypsignathus 
monstrosus, Epomops franqueti, and Myonycteris torquata, which have 
historically existed in a belt across the whole of West Africa. Ge- 
netic analysis performed on samples from 78 individuals in Sierra 
Leone showed no evidence of additional exposure from animals 
sources, beyond the first case. 3 

Phylogenetic analysis 4 estimates the mutation of this strain oc- 
curring sometime around 2002. Recent analysis of samples collected 
at the Lassa Diagnostic Laboratory in Kenema between 2006-2008 
showed that 8.2% (18/219) of patients received with acute illness, 
having tested negative for malaria, tested positive for EBOV (Zaire- 
76 strain) activity. 5 This data should be qualified by the fact that 
the data is a non-random sample. 6 That is, despite being drawn 
from 219 samples, the data may not be indicative of the general 
population. 

A study of immune responses in epidemic and non epidemic 
regions of Gabon was published in 2010. 7 The study sampled 4,349 
individuals in 220 randomly selected villages, finding a 15.5% av- 
erage of specific anti-ZEBOV IgG immune responses. Spatial varia- 
tion of the resistant fraction correlated strongly with landuse, with 
the highest levels (19.4%) found in forested areas. It was speculated 
that the correlation is the result of conduciveness of the environ- 
ment to bats, and that some non-lethal exposure is transmitted 
from fruit. Specifically the authors note in the discussion: "Thus, 
it is possible that ZEBOV antigenic stimulation or aborted infec- 
tion could occur when villagers handle and eat fruits which have 
been contaminated by bat saliva that may contain infectious virus, 



3 Gire et al., 2014 
* Dudas G, 2014 



5 Schoepp et al, 2014 

6 rather than random, it is individuals 
who sought care at a health clinic, 
showed signs of hemorrhagic fever, 
followed a referral, and then tested 
negative for malaria 

7 Becquart et al., 2010 
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inactivated virus, or simple viral antigens. Brief contact with non- 
infectious viral particles is sufficient to induce the type of specific 
immunity observed in our study." The study also noted a linearly 
increasing immune response with age, indicating that children may 
be particularly vulnerable. The antibodies most commonly found 
in the study were active towards the proteins (VP40 & NP), which 
are multi-functional components critical in reproduction and self- 
assembly of the virus. 



3.2 Travel 

On the available mobility within the region, the road network in 
Guinea is well developed by regional standards, both in terms of 
road condition and overall system implementation. The road net- 
works in Sierra Leone and Liberia suffered from neglect during the 
conflict of the 1990s and lingering set-backs to governance, civil- 
society and investment in the decade following. A number of major 
road infrastructure projects were on-going in Sierra Leone and 
Liberia to improve access to rural areas. The outbreak is currently 
showing the most alarming growth in Liberia, in Lofa County and 
the capital Monrovia. Controlling transmission in a mobile popu- 
lation, located in a remote area, with a high level of cross-border 
trade has proven challenging. 

Public transport is typically provided by modified mini-buses 
carrying up to 21 people in space intended for eight. Fuel subsidies 
make travel accessible to average citizens, with travel across an ur- 
ban area in this manner costing as little as $0.25 USD, and regional 
transport (Freetown-Kailahun) costing around $15 USD. 

In rural areas, inexpensive ($1,250) small displacement Chinese- 
made motorcycles have become increasingly ubiquitous. They can 
carry three to four people, market goods and even livestock. A ten 
kilometer trip costs a single passenger about $1.25 USD one way. 



3.3 Compartmental Models 

Basic compartmental models describing the transmission of com- 
municable diseases as grouped populations were first outlined by 
Kermack-McKendrick 8 . 

In brief, compartmental models attempt to quantify the num- 
ber of individuals in various phases of infection with differential 
equations. 

Legrand's Susceptible-Exposed-Infectious-Hospitalized-Funeral- 
Resistant (SEIHFR) model 9 is a compartmental model developed 
specifically for Ebola. Parameters for the transmission of different 
strains of Ebola virus (EBOV) in different historical situations were 
also developed. Of key interest is the transmission coefficients for 
funerals, community activity and hospitals. Specifically, assuming a 
two day traditional burial period, the contact factor for traditional 
funerals was roughly eight times higher than transmission in the 
community and hospitals combined. And that after effective med- 



8 Kermack and McKendrick. 1991 



1 Legrand et al., 2007 
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ical intervention, the transmission coefficients for hospitals and 
funerals is effectively zero. 

4 Working Hypothesis: Endemic-Shift 



Hypothesis: Ebola has previously existed in a stable en- 
demic condition in West Africa, however due to fun- 
damental AND PERVASIVE CHANGES TO THE NATURE OF THE 
TRANSMISSION NETWORK, NEW OUTBREAKS HAVE THE POTEN- 
TIAL TO BECOME LARGE-SCALE EPIDEMICS. 



It is hypothesized here that the epidemic occurred due to the 
wide availability of cellular communication, inexpensive transporta- 



travel, access to multiple trained rural health-care workers and the 
ability of those individuals to seek medical attention in a local hos- 
pital (super-spreading event) were all crucial elements in the initial 
outbreak. They are also factors which have seen dramatic increases 
in recent years. 

In a compartmental model scale, a low reproduction ratio (margir 
above one) in rural areas may indicate that increased mobility and 
population density may be contributing factors. 

If a pre-existing immune response provides some level of resis- 
tance or has a significant effect on mortality, then changes in land 
use and public behavior may also be contributing factors. Given 
the hard-won evidence on the development of immune responses 
(VP40,NP,VP35,GP) in endemic and non-endemic regions 10 , the 
development 11 ' 12 of effective Ebola vaccines based on exposure to 
virus like particles (VLP)s (inducing VP40,NP,GP responses), and 
the effectiveness and safety profile of those vaccines in non-human 
primates (NHP) 13 , some public health recommendations which 
confuse reported or proven pathways of lethal exposure with likely 
sources of VLPs may be in need of review. 

A pre-existing immune response may either be full immunity 
or a key factor in ultimate outcome. Since burials are often crucial 
in the transmission chain, a marginal decrease in fatality has a 
disproportionate on trasnmission effect. 

Four possible scenarios immediately present themselves as tran- 
sitions to a new endemic stability state: 

1. Improved Treatment, 14 increasing the fraction of resistant popu- 
lation members with a more acceptable causality rate. Unless the 
treatment can be administered quickly or self -mediated , it is not 
likely to be scalable. 

2. Vaccination, 15 with priority given to consenting healthcare work- 
ers and their families, and preference to viable options with 
proven efficacy and good safety profiles in NHPs. 



Becquart et al., 2010 

Warfield and Aman, 2011 
Warfield et al., 2007 

Warfield et al., 2014 



14 specifically antibody serums or 
blood transfusions from recent surviors 
and /or early treatment schemes 

15 Galvani et al., 2014 



tion and improved access to health-care facilities. That is, 
to call relatives in neighboring villages, the ability of rela 




? ability 
■s to 
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3. Epidemic, loss of approximately 40%-60% 16 of the current popu- 
lation to the disease itself. 17 

4. An Ebbing Endemic, a wavelike progression of cases controlled 
by a j6 near the critical level, or an oscillating f> reflecting effec- 
tiveness of control measures and changes in public behavior. 

Proper sampling of immunity (by random sampling, capturing 
both the regional variation and uncertainty) is proposed to deter- 
mine the extent of the disease prevalence in the region and the 
required increase in resistant fraction. From the Gabon study 18 , the 
linear increase to age ten, and the leveling of immune response be- 
yond suggests a waning immunity in steady state. From this, it may 
also be necessary to sample individuals in time. Dual species Geo- 
graphic Information System (GIS) compartmental epidemiological 
models have been proposed 19 to study the dynamics of the disease 
with the host species. The development of vaccines with antigen 
specific groups of VLPs suggests techniques capible of characteriz- 
ing specific antigens should be used, such as Western blot analysis. 
The development of raster population density maps and methods 
to quantify movement are a likely by-product of more sophisticated 
analysis being undertaken. 

Assuming that the connectedness of individuals will continue to 
increase as the region is developed, then some fraction of the region 
will require resistance to ZEBOV as some function /(/3), where f> 
will likely continue to increase as the region develops. 



16 rough estimate assuming a period 
of uncontrolled spread, homogeneous 
mixing, followed by complete social & 
economic collapse 

17 the "do nothing" case, or any ap- 
proach which does reverse the course 
of epidemic and cannot scale for 
practical reasons. 
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75.0% 


5.0 
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Table 1: Required immune fraction 

18 Becquart et al., 2010 

19 Olival and Hayman, 2014 



5 SEIR as Partial Quarantine with Pulse Control 

A simplification of Legrand's model for the current situation would 
be to consider a Susceptible-Exposed-Infected/Funeral-Resistant model, 
where quarantined individuals move to the resistant class but a 
fraction proceeds with an unaltered contact rate. The infective class 
here encompasses both individuals who are infected, as well as 
those who succumb to the disease and are buried in a traditional 
manner. The infection rate approaches a limit as the unquarantined 
infected fraction increases. 

Taking both quarantined and unquarantined cases together re- 
sults in a reproduction ratio above zero but less than the uncon- 
trolled rate, 3$ controlled ~ 0/ -^controlled < -^effective < -^uncontrolled- 

Legrand found a basic reproduction ratio of 2.7 (95% CI 1.9-2.8) 
using a stochastic model. However, fitting a deterministic SEIR 
model to the initial outbreak indicates an ^uncontrolled 

of 1.18, 

with an initial 31 effective °f 0-9 after interventions by MSF. This 
may be attributed to the nature of the transmission network and 
be better explained by recent methods developed for hybrid net- 
work /compartmental analysis. 

Following spread to Sierra Leone and Liberia, the ^effective in- 
creased to around 1.24 overall for several months. After the collapse 
of control measures in Monrovia, the regional overall & e ff ec tive nas 



SEIR with Pulse Control Measures: A Model for the 2014 West Africa Ebola Epidemic 6 



increased to « 1.8 — 1.9. This model estimates the progression of 
the infection rate as a logistic regression, where f> approaches some 
value fi max as a function of the infected population. 

Continued transmission implies some unquarantined fraction. 
Unreported cases are not counted in this model. However, the 
unquarantined contact rate together with the current rate of new 
cases (fSSI) may be used to estimate the size of the unreported 
population at some prior time which would lead to a later rate of 
new cases. While, the true extent of unreported cases may be that 
of compounding series 20 , which is not yet treated in this approach. 

In addition to the above assumptions, changes in public behavior 
are difficult to estimate and quantify. And spatial variation can 
have marked effect on contact rates, e.g. between urban and rural 
environments. The progression of the epidemic at different rates 
different environments presents a complex hazard. 

Using the standard Susceptible-Exposed-Infected-Resistant SEIR 
model, with a separate class for deaths, 



1 E~= 0 m 



dS 

It 
dE 
~dt 
dl 

di 
dR 

It 
dD 

~dt 



S': 
E' 
I' ~- 

R' 



-f>Sl, 
PSI - eE, 

al, 
(l-f)al, 



eE 



— D' — foci, 



(la) 
(lb) 
(lc) 
(Id) 
(le) 



S' is change in susceptible 
£' is change in exposed 
where: I' is change in infected 
R' is change in resistant 
D' is change in deceased 

and, 



S is susceptible population, 
E is exposed population, 
I is infected population, 
R is resistant population, 
D is deceased population. 



e is the mean exposed period, expressed in inverse days. 
a is the mean infective period, expressed in inverse days. 
/ is the fatality rate. 

f> is the infection rate proportionality constant, which quantifies 
new infectious contacts per infected person per day. With partial 
patient quarantine and tracing measures in place, the effective 
infection rate is assumed to approach the uncontrolled rate as some 
function of I. A binary operator is added for national quarantine 
measures. 



f> — (\ — Th) x - ^ max , 
p v ; l+me nI 

where: 

T is a binary operator for national house-level quarantine. 
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h is effect of national house-level quarantine. 

f>max is the contract factor with public awareness alone. 

m & n are constants. 

and I (as in indigo) is the infected population. 



A simpler or more accurate estimate of the contact rate can cer- 
tainly be substituted for that component. Time driven estimates 
may be fruitful, however the main driving component appears to be 
the number of people infected. 

The desired effect of house-level quarantine and increased con- 
tact tracing is shown in Figure 2. 




Figure 1: Hypothetical effect of inter- 
mittent house-level quarantine on the 
infection rate and the basic re-production 
ratio, A is the natural frequency of the 
system. 



Time (Qms) 



A substantial increase in contract tracing should cause a dra- 
matic increase in reported cases, this increase should not to be 
confused with an instantaneous increase in infections (from in- 
dividuals exposed five days prior). The converse should also be 
true. That is, if increased contact-tracing reports fewer cases than 
expected, it would provide evidence that enhancing surveillance 
reduces effectiveness. 

Lockdown measures may increase the contact rate on non- 
quarantine days, due to a higher level of economic activity. The 
system has a natural frequency equal to the exposed period plus 
the infective period. Periods of activity and quarantine should be 
staggered across the natural frequency to avoid erratic behavior in 
the underlying system. A weekly schedule is shown on the timeline 
of Figure 2. 

Solving for h may show the effect of house-level quarantine 
on the contact rate. If the results are encouraging, a quarantine 
schedule may be developed to minimize the ultimate size of the 
epidemic, as well as the duration. Even if & e ff ec ti ve does not fall be- 
low 1, mitigating the rate of infection, as soon as possible, increases 
time for development of treatments and vaccines. 

A reliable estimate of p max should be developed to evaluate the 
time-scale of risks posed by 'wait-for-help' strategies. 
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6 Initial Results 

Assuming an infective period of 5.3 days, and an exposed period 
of 5.61 days, the differential equations were fit to total recovered 
individuals and total reported cases for Sierra Leone using sum of 
squared error analysis. 



parameter value 



fimax (total) 
fimax (unit) 

m 
n 



1.81 

4.8e- 

0.29 

0.95 

0.011 



The regression fit is quite sensitive. Particularly, the p max is the 
limit of the curve which can only be expected to be approached, 
but not reached. Using estimate cases may lead to a horrific results. 
As reporting mechanisms approach capacity, there is likely to be a 
breakdown in reporting accuracy. Estimates of 8&t ret from literature 
indicate the rate may in fact be much higher, and even extraordi- 
nary measures may be too late. 



Table 2: Optimized parameters quanti- 
fying breakdown of control measures. 
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Figure 2: Fit on case and recovery data 
to Sept 13. 



The assumed effect of combined national reporting, awareness 
and house level quarantine at given intervals are shown in Figures 
3, 4 & 5. The effect relies on being paired with conventional con- 
tact intervention measures, such as those which have been proven 
successful in more undeveloped countries. 

The true effect can be estimated roughly two weeks after mea- 
sures are taken. Quarantine measures taken in reaction to super- 
spreading events (caused by a lack of security) and implemented 
without contact intervention do not provide a reliable source of 
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data for the effectiveness of the proposed course of action, although 
they may furnish some insight into the uncontrolled rate of spread. 

Figure 3: Effect of 50% effective na- 
Near-term Sierra Leone tional quarantine regiment. Dashed 

lines show 1,2.. .6 days of quarantine 
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Figure 4: Effect of 50% effective na- 
tional quarantine regiment. Dashed 
lines show 1,2. ..6 days of quarantine 
per week. 
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All models are inevitably wrong, the aim is to be less wrong, fur- 
nish useful insight and develop an effective tool for policy makers. 
What policies are implemented, and the manner of their imple- 
mentation, lies with the individuals on the ground those that lead 
them. 
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